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ABSTRACT 

HR8799 is a nearby A-type star with a debris disk and three planetary candidates recently imaged directly. We undertake a coherent 
analysis of various portions of observational data on all known components of the system, including the central star, imaged compan- 
ions, and dust. The goal is to elucidate the architecture and evolutionary status of the system. We try to further constrain the age and 
orientation of the system, orbits and masses of the companions, as well as the location of dust. From the high luminosity of debris 
dust and dynamical constraints, we argue for a rather young system's age of S 50 Myr The system must be seen nearly, but not 
exactly, pole-on. Our analysis of the stellar rotational velocity yields an inclination of 13 - 30°, whereas / S 20° is needed for the 
system to be dynamically stable, which suggests a probable inclination range of 20 - 30° . The spectral energy distribution, including 
the Spitzer/IRS spectrum in the mid-infrared as well as IRAS, ISO, JCMT and IRAM observations, is naturally reproduced with two 
dust rings associated with two planetesimal belts. The inner "asteroid belt" is located at ~ 10 AU inside the orbit of the innermost 
companion and a "Kuiper belt" at i, 100 AU is just exterior to the orbit of the outermost companion. The dust masses in the inner 
and outer ring are estimated to be » 1 x 10"' and 4 x 10"^ Earth masses, respectively. We show that all three planetary candidates 
may be stable in the mass range suggested in the discovery paper by Marois et al. 2008 (between 5 and 13 Jupiter masses), but only 
for some of all possible orientations. For (M(,,Mc,M,() = (5,7,7) Jupiter masses, an inclination / i, 20° is required and the line of 
nodes of the system's symmetry plane on the sky must lie within 0° to 50° from north eastward. For higher masses {Mi,, M^, Mj) from 
(7, 10, 10) to (11, 13, 13), the constraints on both angles are even more stringent. Stable orbits imply a double (4:2:1) mean-motion 
resonance between all three companions. We finally show that in the cases where the companions themselves are orbitally stable, the 
dust-producing planetesimal belts are also stable against planetary perturbations. 

Key words, planetary systems: formation - circumstellar matter - celestial mechanics - stars: individual: HR 8799. 



1. Introduction 

, HR 8799 is an A5V star located at x 40 pc away from Earth, 
' aroun d which three plan etary candidate^ have recently been im- 
' aged (Marois et al.ll2008 ). All three objects have been shown to 
be co-moving with the star. For two of them, also a differential 
proper motion consistent with the orbital motion of a compan- 
, ion around the star has been detected. The prese nce of the outer- 
■ most companion has recently been confirmed bv lLafreniere et al] 
' (|20P9) through analys is of archiva l HST/NICMOS data from 
'. 1998 and by Fukagawa el alj (|2009() with SUBARU/CIAO data 
' from 2002. No further companions with masses greater than 
3Mjup ( Jupiter masses) within 600 AU from the star have been 
found (IClose & MalesI |2009|) . Three imaged companions are 
located at projected distances from 24 to 68 AU from the 
star and presumably have masses of the order of 7-10 Jupiter 
masses, although the mass determination is based solely on 
untested evolutionary models and an assumption that the sys- 
tem's age hes in the range 30-160 Myr. First dynamical analyses 
dFabrvckv & Murrav-Clavl 120081: iGozdziewski & Migaszewskil 
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' For the sake of brevity, we often call them "planets" throughout 
this paper This comes with the warning that the often used definition 
of an "extrasolar planet" as a star-orbiting body with mass under the 
deuterium burning limit remains controversial. Furthermore, it is not 
possible at present to completely exclude the possibility that the mass 
of at least one companion in the HR8799 system lies above that limit, 
although this appears rather unlikely. 



l2009h show that, with the current estimates of planetary masses 
and the assumption that true orbital separations of the planets 
are close to projected distances, the stability of the system on 
timescales comparable to the stellar age is not obvious. Still, the 
system can be stable, as it likely is, for instance if the planets are 
locked in resonances and/or the planetary masses are lower than 
estimated. 

Apart from planets, HR 8799 has long been known to har- 
bor cold circumstellar dust respons ible for excess emission in 
the far-infrared discovered by IRAS dSadakane & Nishidall 19861: 
IZuckerman & SonglllOOl iRhee et alJl2007h . The rather strong 
infrared excess has been a l so con firmed with ISO/ISOPHOT 
measurements dMoor et al.l l2006h . Additionally, Spitzer/IRS 
measurements revealed warrn dust emission in the mid-infrared 
(lJuraetal .1120041: IChen et al.ll2006l) . Both cold and warm dust 
emission must be indicative of one or more dust-producing plan- 
etesimal belts, similar to the Edgeworth-Kuiper belt and possi- 
bly, the asteroid belt in the solar system. Altogether, there ap- 
pears to be an emerging view of a complex, multi-component 
planetary system with several planets, planetesimal belts, and 
dust. 

Given the fact that the planets have been discovered very re- 
cently, it is not surprising that our knowledge of the system is 
very poor and that even the key parameters of the system and its 
components remain vaguely known. A large uncertainty in the 
system's age amplifies the difficulty of inferring accurate masses 
of the companions from evolutionary models, and the mass es- 
timates vary from one model to another even for the same age. 
While there are clear indications that the system is seen nearly 
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pole-on, the exact orientation of its symmetry plane is not known 
either, which makes it impossible to convert projected astrocen- 
tric distances of planets into their true positions. The differential 
proper motion was measured with a reasonable accuracy only 
for the outer planets, and even in that case the accuracy is not yet 
sufficient to constrain orbital eccentricities. As far as the dust is 
concerned, the debris disk remains unresolved, offering no possi- 
bility to better recover the orientation of the system's plane. Even 
the photometry data are scarce. This results in a poor knowledge 
of the dust distribution. 

Obviously, most of these issues could only be resolved with 
new observational data, some of which may come promptly, 
while some others may require longer time spans. However, al- 
ready at the current stage, we find it reasonable to reanalyze 
the av ailable data. While the discovery paper by IMarois et alj 
( I2OO8I) mostly concentrates on th e planets themselves and the 
first "follow-up" p ublications (Fabrv ckv & Murrav-Clavll2008t 
iGozdziewski & M i gaszews ki.2009) provide an in-depth analysis 
of dynamical stability issues, in this paper we try to paint a more 
synthetic view of the planetary system around HR 8799 with all 
its components — central star, planets, and dust-producing plan- 
etesimal belts. 

Section 2 focuses on the stellar properties, notably the stellar 
age (section 2.1) and inclination of the rotation axis, assumed to 
coincide with the inclination of the whole system (section 2.2). 
In Section 3, we analyze the observed spectral energy distribu- 
tion (SED, section 3.1) and try to fit the data with several dust 
belts (section 3.2). Section 4 treats presumed planets and tries to 
constrain their masses both from evolutionary models (section 
4.1) and dynamical stability requirement (section 4.2). Section 
5 checks whether planetesimal belts, as found to fit the infrared 
photometry, would be dynamically stable against planetary per- 
turbations. Section 6 contains our conclusions and a discussion. 



2. The central star 

2.1. Age 

iMarois et al.1 (|2008|) give an age of 30 to 160 Myr for HR 8799 
and hence its companions, considering several age indicators. 
While most of the indicators support a rather younger age, few 
of them still allow for an older age or even suggest it. One indica- 
tor is the galactic space motion (UVW) of the primary, which is 
similar to that of c lose young associ ations (Marois et al. 2008}). 
Using these data, iMooretalJ (|2006|) propose HR 8799 to be a 
member of the Local Association at an age of 20 to 150 Myr 
with a probability of 62 %. Another method is the position of 
HR 8799 in the Hertzsprung-Russell diagram. Taking into ac- 
count the low luminosity of the star (after co rrection for its low 
metallicity) as well as the UVW space motion. lRhee et aL I (f2007h 
arr ived at an a ge of 30 Myr 

IMarois eF al. (2008) further note that A Bootis stars are gen- 
erally thought to be young, up to a few 100 Myr. However, the 
Hipparcos mission has shown that the well established A Bootis 
stars of the Galactic field comprise the whole range from the 
Zero Age Main Sequence to the Terminal Age Main Sequence, 
which is ~ 1 Gyr for an A-type star (Turcotte 2002; Paunzen 
l200li and references therein). The best indicator for a rather 
older age is the location of HR 8799 in a T^s vs. \ogg diagram 
derived from pub lished uvby/J photometry. Using this method, 
ISong et all (1200 ll) find an age of 5 to 1128 Myr with a best 
estimate of 732 Myr and lChen etlD (2006) an age of 590 Myr. 

An independent argument in favor of a rather younger age 
may come from the the dust portion of the system. The mea- 



sured infrared excess ratio of ~ 100 at 60 -90 //m (see Fig. [T] 
below) would be typical of a debris disk star with the age of 
$ 50 Myr (see Su et al. 2006, their Fig. 5). However, this ar- 
gument is purely statistical and must be taken with caution. For 
instance, one cannot exclude the possibility that the formation of 
such a planetary system with three massive planets in very wide 
orbits could trace back to an exceptionally dense and large pro- 
toplanetary disk. The latter might leave, as a by-product, a more 
massive debris disk on the periphery, whose fractional luminos- 
ity might be well above the statistically expected level. 

Altogether, there seem to be more arguments to advocate 
a younger age of the system of the order of several tens of 
Myr On any account, as pointed out in the discovery paper by 
IMarois et al.l ( l2008l) and is further discussed in section 4 below, 
extremely old ages would inevitably imply high object masses in 
the brown dwarf ra nge — for all the evolutio r iary models used to 
infer the masses. Fabr vckv & Murrav-Clavl (1200 8) have shown 
that dynamical stability of such a system is problematic. It can 
be stable for masses up to 20 Jupiter masses, but only for very 
special orbital configurations. 

2.2. Rotational period and inclination 



As summarized by ISadakanj (|2006|) the Vega-like, y Doradus 
type pulsator HR 8799 is showing A Bootis-like abundances. He 
concludes that for the case of HR 8799, which is known to be a 
single star and associated with a dusty disk, the scenario invok- 
ing the process of selective accretion of circumstellar or inter- 
stellar material depleted in refractory elements, seems to be the 
favorable explanation for the unusually low abundances. 

The fact that HR 8799 is a y Doradus type pulsator makes 
it hard to find an indisputable rotational pe riod. However, sev- 
eral authors give a value of ~ 0.51 days (e.g. 'Rodriguez & Zerbil 
1995). Zerbi et al. ( 1999) found from a multisite campaign three 
independent frequencies (0.5053 d, 0.5791 d, 0.6061 d) and a 
coupling term between them (4.0339 d). All of these frequen- 
cies could be independent pulsational modes. However, if one of 
these frequencies corresponds to the rotational period of the star, 
we are able to calculate the inclination / of its rotational axis. 

From the possible rotational frequencies (0.5053 d, 
0.5791 d, 0.6 061 d) and from the r a dius of HR 8799 
of 1 .32 j?B (lAllende Prieto & Lamber g Il999l) to 1.6 Rq 



(iPasinetti Fracassini et al.l 1200 ll) . we determine the possi- 
ble range of the true rotational velocity v of the star of 
1 10-160 kms '. These values agree quite well with the median 
vsin; of A4— A6 main sequence stars of 159 + 7.2 km s"' 
(iRover et al.ll2007l) . On the other hand, the projected rotational 
velocity v sin / of HR 8799 was measured by several authors to 
be be t ween 35.5 km s~' and 55 km s~' (e.g. lKaye & Strassmeied 
Il998t IUesugi&Fuk"udalll982l) . From V and vsin/, we finally 
derive a possible range of the inclination of the star of 13°-30°. 
Note that in the above estimates we excluded the 4.0339-day 
period as this would lead to sin / > 1 . 

It can be expected that the rotational equator of the star 
and the planetary orbits are all aligned with each other. Spin- 
orbit alignment is a common assumption, consistent for instance 
with the data of most transiting planets (e.g. Fabrvckv & Win^ 
200 91). It has also been confirmed for Fomalhaut and its disk 
(iLe Bouquin et al.l 120091) . Nonetheless, a misalignment on the 
order of several degrees is likely. It is exemplified by our own 
system, in which the Jupiter orbit is tilted by ~ 7° to the solar 
equator and the orbital planes of giant planets differ from each 
other by a few degrees. With this caveat, we make an assumption 
of the perfect alignment throughout the paper This assumption 
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is particularly needed to reduce the complexity of our dynamical 
stability studies presented below. 

An additional argument in favor of a low inclination of the 
whole system comes from the fact that the measured proper mo- 
tion of the companions relative to the star is also consistent with 
a nearly pole-on view onto the system (Marois et al. 2008). The 
most recent astrometric analysis by Lafre niere et al.l (1200 9) fur- 
ther supports this conclusion. Their best fits with circular or- 
bits suggest the semimajor axis a ~ 68-74 AU and the incli- 
nation / ~ 13°-23° for the orbit of the outermost companion, 
HR 8799 b. 

3. Dust and planetesimal belts 

3.1. Observed SED 

Table[T]lists the catalogs and references used to provide the opti- 
cal, infrared, sub-mm and millimeter photometry. We employed 
the Hipparcos and Tycho databases as well as USNO and GSC 
catalogs to compile the optical photometry, whereas the near- 
infrared data were taken from the 2MASS survey. The mid- and 
far-infrared photometry is provided by the IRAS and ISO satel- 
lites, while sub-mm and millimeter data were obtained at JCMT 
and IRAM. For transforming the B, V, R, I magnitudes into units 
of flux [Jy] we used the standard calibration system of Johnson, 
whereas in case of t he 2MASS J, H, Kg bands the calibrations of 
ICohen et al.l (12003 ') were applied. For t he IRAS fluxes a color- 
correction factor ^eichmanet aL|ll988|) assuming a black body 
spectral energy distribution for a temperature of 5000 K was em- 
ployecfl 

The optical and near-infrared photometry was used to de- 
rive a best-fit photospheric model. We estimated the stellar pho- 
tospheric fl uxes by minimum fitting of NextGen model at- 
mospheres dHauschildt et al.iri999l) . using only bandpasses with 
wavelengths shorter than 3 fim, where no excess emission is ex- 
pected. In our search for the best fit we employed a system of 
NextGen models with an effective temperature step size of 200 
K, with a log g of 4.5, and with solar metallicit}0 Varying the 
temperature as well as the stellar radius, which aff'ects the solid 
angle dilution factor, we derive a best-fit temperature of 7400 K 
and a best-fit radi us of 1.34 Rir,, both i n a very good agreement 
with the results of lGrav & Kav3 (Il999l) . 

In order to obtain more information on the dust compo- 
nent of the system, we extracted publicly available IRS data 
for HR8799 from the Spitzer archive with the Leopard soft- 
ware ( http://ssc.spitzer.caltech.edu/irs/). Those data were taken 
in Decem ber 2003 under AOR 3565568 and originally published 
bv lChen et al., t2006t) . Using the post-basic correction data pro- 
vided by the standard IRS pipeline, we first subtracted the zo- 
diacal light background estimated by the SPOT software. Then 
we joined the datasets for the individual IRS modules and orders 
and subtracted a simple Rayleigh- Jeans stellar photosphere in 
the wavelength region of interest, assuming now that there is no 
excess at 8 i-im where the IRS spectrum starts. For plotting and 
fitting purposes, the scatter in the resulting dust spectrum was 
reduced by taking averages over 4 neighboring data points. This 
increases the signal-to-noise ratio without lowering the spectral 
resolution, which remains constrained by the instrumental value 
R ~ 60-120. Spurious features at the edges of the spectral orders 
around 14 and 20 yum were removed. 

^ Correction factors are only available over a coarse temperature grid 
. . . , 5000 K, 10000 K. 

3 Although the metallicity of HR8799 is -0.47 jGrav & Kavell999l) , 
this parameter has no important effect on the fitting result. 
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The resulting spectrum shows a somewh at (by ^ 40%) 
weaker excess than that by IChen et aTl ( l2006l) . especially be- 
tween 20 and 30 fim. One possible origin of this discrepancy 
is the approach used for background subtraction. In addition, 
both spectra fall well below the IRAS upper limits at 25 yum 
but also by several cr below the IRAS 12 fim faint-source and 
point-source measurements. The easiest way to reduce this scal- 
ing uncertainty would be using a good photometry point in the 
mid-infrared. Unfortunately, a 24 //m Spitzer/MIPS point, which 
would ideally serve for this purpose, is not yet available in the 
Spitzer archive, although these data have recently been taken 
(Kate Su, pers. comm). Whether the maximum at around 1 1 jum 
is a silicate feature is unclear. Definite conclusions would require 
a longer exposure time and a higher spectral resolution. 

The resulting set of photometry points and the IRS spectrum 
are shown in Fig. [T] (with photosphere) and in Fig. |2] (excess 
emission only). 

Table 1. Photometry of HR8799. IRAS fluxes are color cor- 
rected as described in the IRAS Explanatory Supplement 
dBeichman et al.lll988l) . 



Photometric band 


Flux or magnitude 


F,„„, (IRAS) 


Ref. 


[mag] 


B 


6.090 ± 0.300 




(1) 


B 


6.196 




(2) 


B 


6.210 ±0.010 




(3) 


B 


6.214 ± 0.009 




(4) 


V 


5.960 ± 0.010 




(4) 


V 


5.959 




(2) 


V 


5.960 ± 0.010 




(3) 


R 


5.810 ±0.300 




(1) 


I 


5.690 ± 0.300 




(1) 


J 


5.383 ± 0.027 




(5) 


H 


5.280 ±0.018 




(5) 


K, 


5.240 ±0.018 




(5) 


[Jy] 


IRAS PSC 12 iim 


0.267 ± 0.034 


3 


(6) 


IRAS PSC 25 iim 


0.246 ± 0.000 


1 


(6) 


IRAS PSC 60 ;im 


0.307 ± 0.061 


2 


(6) 


IRAS PSC 100 fim 


2.376 ± 0.000 


1 


(6) 


IRAS FSC 12 fim 


0.278 ± 0.036 


3 


(7) 


IRAS FSC 25 fim 


0.174 ±0.075 


1 


(7) 


IRAS FSC 60 fim 


0.311 ±0.062 


3 


(7) 


IRAS FSC 100 fim 


3.202 ± 0.977 


1 


(7) 


ISO 60 fim 


0.412 ±0.021 




(8) 


ISO 90 fim 


0.585 ± 0.041 




(8) 


[mJy] 


JCMT 850 fim 


10.3 ± 1.8 




(9) 


JCMT 1 100 Mm 


< 33 




(10) 


IRAM 1200 


4.8 ± 2.7 




(10) 



Referen ces: (1) TheUSNO-Bl .O Catalog 4Monet et aU2003l); (2) NOM AD 
Catalog iZ acharias et alj2004ft . from Tycho-2 Catalog lH0s et al. 2000); (3) 
The Guide Star Catalog, Version 2.3.2 07asker et al. 2008); (4) The Hippai'cos 
and Tycho Catalogues I Perrvmair& ESA 1997); (5) 2MASS All-Sky Catalog 
of Point Sources (Skrutskie et al. 2006); (6) IRAS catalogue of Point Sources, 
Version 2.0 iHelou & Walker 1988); (7) IRAS Faint Source Catalog, \h\ > 10, 
Version 2.0 i Moshir et al. 1990); (8) (Moor et al. 2006) (9) 
IWilliams & Andrewsi200a) (10) ISvlvester et all 19961) 



3.2. Interpretation 

To get a rough idea of the location of the dust belt(s) and the 
amount of dust in the HR8799 system, we modeled the SED 
assuming a double power-law surface number density of dust 
oc s^'^r^^, where s and r are the dust grain size and distance 
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1 10 100 1000 

wavelength [nm] 

Fig. 1. Spectral energy distribution of HR 8799. Circles are pho- 
tometric measurements, triangles are upper limits (TablelTJ. The 
black line between w 5 and 40 fim is the Spitzer/IRS spectrum. 
Overlaid is the best-fit NextGen model (grey line) for the stel- 
lar photosphere. Its estimated relative error is plotted under the 
main panel. 

Table 2. Locations and names of the first-guess dust rings. 



ring location 


ring extension [AU] 


name 


inside d 


3-15 


ring d 


between d and c 


28-32 


ring cd 


between c and b 


45-60 


ring be 


outside b 


75 - 125 


ring b 



from the star, respectively. Keeping in mind that SED interpre- 
tation is a degenerate problem and to decrease the number of 
free parameters, we restrict ourselves to the case of ^ = 3.5 and 
^ = 1. The dust composition was assumed to be astronomical 
silicate (Laor & D rainel[T99 3). As the minimum grain size, we 
chose the radiation pressure blowout radius, which amounts to 
X! 5 fim for astrosilicate. The maximum grain radius, which has 
little effect on the results, was arbitrarily set to 1000 /jm. 

At first, we modeled the SEDs that would arise from four hy- 
pothetical dust rings with arbitrarily chosen extensions, located 
inside the orbit of the innermost planet HR8799 d, between d 
and c, between c and b, and outside the orbit of the outermost 
planet HR8799 b (Table |2]l. Dust masses were chosen in such a 
way as to reproduce the 60 //m flux. The results are shown in 
Fig.|2] 

Comparison of the first-guess model SEDs with the available 
photometry and spectrometry observations reveals two issues. 
First, one single ring is not capable of reproducing the entire 
set of observations from the mid- to the far-infrared. Second, if 
the 10 jum silicate feature in the IRS spectrum is real, a substan- 
tial fraction of particles smaller than the blowout size will be 
required. 

Taking these discrepancies into account, we combined two 
rings, 'ring d' and 'ring b', and fitted the lower size cutoffs and 
dust masses to the observations. In our 'best fit' model, which is 
shown as solid line in Fig. |2] the minimum sizes are 2 and 6 fim 
and the dust masses are 1.4 X 10"^ Me and 4.2 x 10"^ foj. 
the inner and outer ring, respectively. As a word of caution, we 
remind that these estimates directly inherit the uncertainty of the 



IRS spectrum calibration, as discussed in section 3.1, which is a 
factor of several in dust mass. 

A question arises whether these fits are physically rea- 
sonable. In particular, it is not immediately clear whether the 
blowout particles are really needed to reproduce the observa- 
tions. Similar problems have already been enco untered in stud - 
ies of other debris disks (e.g. around Vega, see ISu et alJl2005h . 
However, such inconsistencies can be overcome by consider- 
ing a complete dynamical treatment of the specific disk with- 
out significant changes in the disk location (Miiller et al., in 
prep.). Further, our choice of power-law spatial and size dis- 
tributions is only a rough, albeit commonl y used, approxima- 
tion. More elaborate dynam ical studies (e.g.l Krivov et al.|[2006l : 
iThebault & Augereaull200'7l) show clear deviations from this as- 
sumption, especially in the case of the dust size distribution 
where a wavy pattern arises from an underabundance of small 
particles induced by radiation blowout. Thus, we checked the 
impact of moderate variations in the slopes (0.1 < ^ < 1.9 and 
2.5 < q < 4.5). It turned out that the slope of the spatial dis- 
tribution ^ has little effect on the results, moderately changing 
only the short-wavelength part in the SED from the inner dust 
ring. The size distribution slope q, however, affects the resulting 
emission appreciably. While for the inner ring the changes are 
still small (a steeper slope would amplify the silicate features 
at 10 and 20jum whereas a shallower distribution would wash 
them out, just as expected), similar changes in the outer ring 
would require a strong compensation by altering other disk pa- 
rameters (which were fixed in our approach). The reason is that 
for the outer component both the rise and the fall of the SED are 
well constrained by photometric observations. They place tight 
constraints on the width of the SED, in contrast to the inner ring 
where only the short- wavelength part of the SED is known. Since 
dust in the outer ring is much colder than in the inner ring it is 
no longer the strength of the features but the width of the SED 
that is affected by a different slope in the size distribution: the 
steeper the distribution, the narrower the SED. On any account, 
for our purpose of determining the rough location and mass of 
dust the simple fitting approach used here is sufficient. 

Another question is, how well the edges of the outer and in- 
ner ring are constrained. To check this, we varied them and fitted 
the SED again, leaving the lower cutoff size and dust mass as free 
parameters. For the outer ring, we found reasonable fits with the 
inner edge between 75 and 120 AU and the outer edge between 
125 and 170 AU (for q - 3.5 and ^ = 1). The outer edge in the 
inner ring can go from 15 down to 10 AU. In fact, the inner ring 
truncated at 10 AU provides even slightly better agreement with 
the IRS spectrum between 20 and 30 ;um. However, due to the 
calibration uncertainties it is difficult to assess the accuracy of 
the fit, which leaves the outer edge of the inner ring rather un- 
constrained. The inner edge of the inner ring can be as close to 
the star as 2 AU to conform to the IRS spectrum. 

So far we have discussed dust rings in the system. However, 
the presence of a dust belt requires a belt of planetesimals that do 
not show up in the observations, but produce and sustain visible 
dust. Due to the radiation pressure effect, dust grains typically 
move in orbits with periastra at the planetesimal belt and apas- 
tra outside it. The smaller the grains, the farther out from the 
star they are spread. Thus the dust-producing planetesimal belt 
is expected to be narrower than the observed dus t ring and to 
lie within the dust ring close to its inner edge (e.g. lKrivov et al.l 
i2006l) . Therefore, we can expect that the outer planetesimal belt 
is located at ~ 75-120 AU and the inner planetesimal belt at 
2-10 AU from the star It is important to check whether the ex- 
pected locations of the outer and inner planetesimal belt are dy- 
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ring d 
ring cd 
ring be 
ring b 
best fit 




dynamical evolution of HR8799 with 1=0° 



200 



Fig. 2. Excess emission of HR 8799 in the infrared. The grey 
sohd line is the IRS spectrum. The circles are the observed 
fluxes, and the triangles are upper limits. Dotted, short-dashed, 
dash-dotted, and long-dashed lines are the SEDs from the four 
first-guess dust rings (Table |2]i. The black solid line is our 'best 
fit'. 



namically compatible with the presence of the outermost and in- 
nermost planet, respectively. We will investigate this in section 
5, trying to find additional constraints on the location of the two 
rings. 



4. Planets 

4.1. Masses from models 

A common method to estimate masses of various astrophysical 
objects, from stars to planets, is to use their formation and evolu- 
tion models. Such models predict essential physical parameters 
of objects of various masses, notably their luminosity and tem- 
perature, as a function of age. If the age is known, a comparison 
with luminosity or temperature retrieved from observations al- 
lows one to estimate the mass. In the case of HR8799, the tem- 
perature of the companions is unknown, but their luminosity has 
been derived from brightness and distance with suflicient accu- 
racy to apply the models. Ages, masses, and other parameters of 
ot her directly imaged p lanet candidates were recently reviewed 
bv lSchmidt etaTI J2009l) . 

However, all models involve simplifying assumptions and 
adopt certain initial conditions, e.g. the initial internal energy 
and temperature structure. Further, the so-called hot-start mod- 
els all start at a non-zero age (e.g., 0.1 or 1 Myr) with a finite 
luminosity and thus do not cover the actual formation stage. For 
these reasons, models may n ot deliver reUable results for at least 
the first several Myrs (see Wuchterl 2001; Chabrier et al. 2005, 
for discussion), if not several hundred Myrs ( Stevenson 1982 ). 

Using the model by iBaraff'e et al.1 (l2003h an d assuming the 
age range of 30-160 Myr, iMarois et al] (l2008h estimated the 
masses of companions to be 7Mjup (5-llMjup) for HR8799 b 
and lOMjup (7-13Mjup) for HR 8799 c & d. We have recalculated 
possible masses of the three companions with the aid of several 
state-of-the-art hot-start evolutionary models from system's age 
and companions' luminosity. In so doing, we allowed a broader 
range of possible ages, as discussed in section 2.1. Table[3]shows 
the results. Note that some models give masses onl y for limited 
age and/or mass ranges; for instance, the model bv iMarlev et alj 




100 
time [kyr] 



200 



Fig. 3. Instability of the planetary system with nominal masses 
in the non-inclined case. For each of the planets, three curves 
correspond to semimajor axes, pericentric distance, and apocen- 
tric distance. 



(I2007h does not go beyond lOMjup. This explains why some po- 
sitions in the table are not filled. 

On t he whole, the mass estimates we obtain are similar to 
those bv lMarois et al. I (120081) . The ages above ~ 160 Myr would 
lead to "non-planetary" masses above the standard deuterium 
burning limit of 13.6Mjiip. On the other hand, none of the models 
predict masses below 3-5 (HR8799 b) and 6 - 8 (HR8799 c, d) 
Jupiter masses even for the youngest plausible ages of ~ 20 Myr. 

4.2. Masses from stability requirement 

The simplest possible assumption one can make is that the sys- 
tem is seen perfectly face-on (inclination / = 0) and that all 
th ree planets are initially in circular , co-planar orbits. However, 
as Fabrv ckv & Murrav- Clav (2008) pointed out, such a system 
with masses reported in Marois et al. (2008) would be unsta- 
ble. This is readily confirmed by our numerical integrations 
that are described below. Fig. [3] shows the time evolution of 
the planetary orbits; already after 134 kyr a close encounter be- 
tween HR8799 c and d occurs and triggers instability of all 
three planets ending with the ejection of HR8799 c 50 kyr later 
These timescales are fully con sistent with those reported by 

iFabrvckv & Murrav-Cla v* ('2008V 

The studies by [F abrvck v & Murrav-Clavl (l2008l) and 
iGozdziew ski & M igaszewskil 020091) were largelv based on fit- 
ting simultaneously the observed positions and differential 
proper motion of the companions and checking the stability of 
the resulting systems in the course of their dynamical evolu- 
tion. Noting that constraints on eccentricities from the differen- 
tial proper motion are rather weak, but clear indications exist for 
inclined orbits from the rotational period analysis (see section 
12.21 ), here we employ a different method. We confine our simula- 
tions to initially circular orbits (and let the eccentricities evolve 
to non-zero values at later times), but allow the symmetry plane 
to have all conceivable non-zero inclinations and an arbitrary 
orientation. Thus, for the subsequent analysis we introduce two 
angles (Fig.|4|. One is the inclination / itself, measured between 
the angular momentum vector and a vector pointing toward the 
observer. Another angle is the longitude of node Q of the sys- 
tem's symmetry plane on the plane of the sky, which is measured 
from north in the eastern direction. We vary the inclination / from 
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Table 3. Masses of HR 8799 b / c / d from luminosity (and absolute K-band magnitude) using various evolutionary models. 



HR8799 b 


Luminositv los L/L„ = -5.1 ±0.1 (Marois et al. 2008) 


Model 


20 Myr 


30 Myr 


60 Myr 


Mass [Mj„p] at age 
100 Myr 160 Myr 590 Myr 730 Myr 


1000 Myr 


1128 Myr 


Burrows et al. (1991) 


3.5-4.5 


4.5-6 


7-8.5 


9- 11 11.5- 12.5 22-26 25-30 


28-33 


30-36 


IVlallCy CL ai. \jl,\J\J/ ) 


3-5 


4-7 


6-10 








cnaurier et ai. (zvuv) 








21-26 


30-35 




oaraire et ai. [zxjuj) 




4-5 


6-7 


9-10 21-26 


30-35 




Baratte et al. (zUUs) 






- 7 


-9 






Baratte et al. (zUUi) 




-5.5 


- 8.5 


- 10.5 - 30 


-38 




HR8799 c / d 


Luminositv log L/L„ = -4.7 ± 0. 1 (Marois et al. 2008) 


Model 


20 Myr 


30 Myr 


60 Myr 


Mass [Mj„p] at age 
100 Myr 160 Myr 590 Myr 730 Myr 


1000 Myr 


1128 Myr 


Burrows et al. (1997) 


6-7.5 


7.5-9.5 


11 - 12 


12.5- 13 13- 13.5 30-38 35-43 


40-48 


41-50 


Marlev et al. (2007)" 


6-8 


8-10 










Chabrier et al. (2000) 




6-7 


8-10 


10-11 28-34 


39-46 




Baraffe et al. (2003) 




6-7 


8-10 


10-11 27-31 


37-43 




Baraffe et al. (2008)* 






-9 








Baratfe et al. (2003)'' 




-7.5 


- 10.5 


-11.5 -38 


-48 





Remarks: All mass values were shortened to appropriate deci mals and partly rou nded to halves or integers of a Jupiter mass, (a) Hot-start models, 
(b) Non irradiated models, (c) Using Mjf ,= 14.05 ± 0.08 mag iMarois et alj2008 ). (d) Using Mk., - 13.1 2 mag instead of slightly dilferent values 
Mjf, = 13.13 ±0.08 mag for HR8799C and Mkj=13.11 +0.12 mag for HR8799 d, given in Marois et alj ( l2008l) . 




Fig. 5. Semimajor axes of HR8799b, c and d (solid lines) and 
their Hill radii for nominal masses (dashed) as a function of Q. 
Fig. 4. Orientation of the system with respect to the line of sight, for / - 30°. 



0° to 45°, thus extending the range suggested by section 2.2 to 
higher values. This may be useful to accommodate a possible 
tilt of planetary orbits to the stellar equator The rotation angle 
Q is unconstrained by the observations. It is sufficient to vary it 
from 0° to 180°, since the true mutual positions of all three plan- 
ets at (i, Q) and (/, Q. + 180°) would be exactly the same. For 
each (/, Q)-pair we can convert the observed (projected) instan- 
taneous positions (astrocentric distances and positional angles) 
of the three planets into their true positions in space. If, further, 
initially circular orbits are assumed, the calculated distances of 
planets will coincide with their initial semimajor axes. We thus 
consider a two-parametric (/, Q) set of possible systems; for each 
of them, the initial orbital configuration is fully and uniquely de- 
fined. 

Fig. |5] shows the astrocentric semimajor axes of the three 
planets with an inclination of / = 30° as a function of rotation 
angle Q., as well as the corresponding Hill spheres for nominal 
planetary masses. Similarly, Fig. |6] depicts the initial semima- 
jor axis of the planets depending on / and Q. Note that the recent 



analysis of HR8799 b by lLafreniere et all (l2009l) . which is based 
on the archival HST/NICMOS data from 1998, yielded a semi- 
major axis of fl ~ 68-74 AU and an inclination of / - 13°-23°. 
Finally, Fig. |7] plots the difference of initial semimajor axes of 
HR8799 b and c, as well as HR8799 c and d, again depending 
on / and O. From all these figures, it is clearly seen that the or- 
bital spacing, and therefore the stability, might indeed strongly 
depend on the orientation of the system. 

This expectation is fully confirmed by the main bulk of 
numerical integrations that we performe d with the aid of the 
MERCURY6 package ("Chambers" 1999"). We used the hybrid 
symplectic/Bulirsch-Stoer integrator with an adaptive stepsize 
and a 10"'"^ angular momentum conservation accuracy, which 
changes to Bulirsch-Stoer algorithm at distances less than 3 Hill 
radii. Output has been stored every 1000 years. Each integra- 
tion terminated when two planets had a distance less than half 
the Hill radius or after an integration time of t^^^ = 100 Myr. 
In all cases we assumed a stellar mass of 1 .5Mq and a distance 
of 39.4pc to convert the separation angle into the projected as- 
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Fig. 7. Difference of initial semimajor axes Aa between 
HR8799 b & c (top) and HR8799 c & d (bottom) as a function 
of / and Q., as well as the location of the 2:1 commensurability 
of the periods for circular orbits. 

Table 4. Three sets of planet masses used in numerical integra- 
tions 



HR8799 


b 


c 


d 


low mass 
nominal mass 
high mass 


5Mj„p 
7Mj„p 
llMj„p 


7Mj„p 
10Mj„p 
13Mj„p 


7Mj„p 
10Mj„p 
13Mj„p 



Fig. 6. Initial semimajor axes of HR8799b (top), c (middle) and 
d (bottom) as a function of / and Q. Solid horizontal lines in the 
uppermost panel border the ran ge of inclination range o f i - 13- 
23° for HR8799 b, reported by iLafreniere et al.1 (|2009|) . Dashed 
lines do the same for the semimajor axis range, a - 68-74 AU. 



trocentric distance. The three planets s tarted at positions a t the 
epoch of 2008 Sept. 18 (see Table 1 in lMarois et ani2008h . We 
used 3 different sets of possible planet masses (see Table |4]i and 
checked stability for / e {0, 1,...,45°) and O € {0, 5, . . . , 180°). 
Note that, although we assumed initially circular orbits, it does 
not mean that the orbits stay circular at later times. Just the oppo- 
site: the mutual perturbations always force eccentricities to take 
values in the range between zero and approximately 0.1, so that 
the initial circularity is "forgotten" by the system. 

The results of these integrations are presented in Fig. |8] It 
depicts the time interval until the first close encounter — as a 
proxy for stability — with an upper limit of lOOMyr It is seen 
that all three planets may be stable for 100 Myr for either set 
of planetary masses, but only for some of all possible orienta- 



tions. Specifically, an inclination of ^ 20° is required and the 
rotation angle Q. must lie within the range from 0° to 50°. 
The higher the masses, the narrower the "stability spot" in the 
(/, f2)-plane. Further, we conclude that it is the inner pair (c and 
d) that tends to destroy the stability. Indeed, comparing Fig.[8]to 
Fig. |2] one sees that the most stable regions are those where Aa 
of HR8799 c and d is the largest, whereas that of HR8799 b and 
c is not. Finally, a comparison of Fig.[8]and Fig.|6]shows that the 
position of the "stability spot" in Fig.[8]roughly matches the in- 
clination of / = 13-23° and the sem imajor axis of a = 68-7 4 AU 
of the outermost planet reported by ILafreniere et all (l2009h . 

We now check whether all considered geometries are con- 
sistent with the measured differential proper motion of the com- 
panions. Fig. |9] depicts the projected differential proper mo- 
tions /i that the planets in circular orbits would have for each 
pair of / and Q.. Overplotted are the values of p actually mea- 
sured together with their Icr and 2cr deviations. For HR8799 c, 
we used p - 30 + 2mas/yr from Marois et al. (2008). In con- 
trast, for HR8799 b we derived the differential proper motion of 
22 ± 2 mas/yr by combining the Marois et al. measurements with 
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Fig. 8. Duration of stability depending on / and Q for different 
sets of planetary masses given in Tab. H) top: low, middle: nom- 
inal, bottom: high. A grid of dots corresponds to the actual set 
of numerical runs (one dot - one run). Curves show where the 
periods between b-c and c-d would have a commensurability of 
2: 1 if the orbits were exactly circular. 



the 1998 data from 'L afreniere e t al.' (lO QSj). N ote that further 
adding the 2002 positions from Fukagawa et al.l (I2009.) does not 
change this result. For co mparison, the differ ential proper mo- 
tion of HR8799 b given in' Marois etli\ (l2008h is 25 ± 2 mas/yr 
In the non-inclined case, the calculated differential proper mo- 
tion is always within Icr of the measured one. For HR8799 b, 
all orbits with an inclination < 33° lie within Icr of the mea- 
sured value. For HR8799 c, the same is true for / < 28°. Nearly 
the whole parameter range of (/, Q) explored here is compatible 
with observations within 2cr. It is easy to show that taking into 
account actual low eccentricities up to » 0.1 acquired by the 
planets would not change this conclusion. 




60 80 100 120 
rotation angle [deg] 

differential proper motion HR8799C (n = 30 ± 2 mas/yr) 
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Fig. 9. Sky-projected differential proper motion that HR8799 b 
(top) and c (bottom) would have for each pair of / and Q. (as- 
suming circular orbits). The isolines are atp,p- cr, and p - 2cr 
{(7-2 mas/yr), where p is the differential proper motion de- 
rived from the measurements, as described in the text. Note that 
the p + cr and p + 2cr curves for both companions fall outside the 
plots. A panel for HR8799 d is not included because the uncer- 
tainty ip - 42 mas/yr +27 mas/yr) is too large. 



To understand the stability region shown in Fig.|8] we check 
the possibility that it may be related to resonances. From Fig.fTOl 
it is obvious that the two inner planets are close to the 2:1 res- 
onance. With still coplanar and circular orbits, we looked for 
combinations of /, f2 that would correspond to the 2:1 commen- 
surability. The resulting loci of the nominal 2:1 resonance in 
the (/,Q)-plane (Fig. [8]l just encircle the stability region. This 
strengthens the hypothesis that the stability may be directly re- 
lated to the 2:1 resonance. 

We then checked whether or not the orbits of HR8799 c 
and d within the stability region are indeed locked in the res- 
onance. To this end, we have calculated the resonant argument 
i^cd = 2/lc - /Id - (jL>i, where Ac and Ad are the mean longitudes 
of HR8799 c and d and is the argument of pericenter of the 
latter planet. We found that all stable orbits are indeed resonant. 
Interestingly, the initial values of Ac and Ad adopted in all nu- 
merical runs were such that the planets are not locked in the 
resonance initially, but — in all stable cases - the system swiftly 
"slips" into the resonance, which makes it safe. The resonant 
argument i/Pcd librates around 0° with an amplitude (which we 
calculated as a standard deviation) of 22°-100°. For compari- 
son, the non-resonant case would have a standard deviation of 
~ 103.9°. However, we made sure that even the cases with the li- 
bration amplitude up to 100° are resonanceces, albeit shallow. In 
these cases, the resonant argument circulates rather than librates. 
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Fig. 10. Semimajor axes of HR8799b, c and d (solid lines) and 
positions of external 2:1 (dashed) and 3:2 nominal resonances 
(dotted) of two inner planets as a function of Q for / = 30°. 
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Fig. 11. Typical behavior of resonant arguments for the 2: 1 mean 
motion resonance between HR8799 b & c {top) and HR8799 c 
& d {bottom). Nominal planetary masses are assumed. The libra- 
tion amplitude is 86° {left top), 36° {left bottom), 63° {right top), 
and 60° {right bottom). 



Thus the phase trajectories on the cos (pd-ec sin (f^d plane are 
circles with an offset from (0,0), which is indicative of a resonant 
locking. 

For all stable configurations we then calculated the resonant 
argument for the two outer planets, HR8799 b and c, defined 
as (|Pbc = 2/lb - Ac - ol>c. We found that they are in a 2:1 res- 
onance, too, with the standard deviation in the range 36°-97°. 
Note that in all stable cases at least one of the resonances, c-d 
or b-d, is strong, as suggested by a low libration amplitude. Two 
typical examples of the time evolution of the resonant argument 
for both planetary pairs are shown in Fig. fTTl Thus our results 
are consistent with those by [pa brvckv & Murrav-Clavl (l2008l) 
and lGozdziewski & Mi gaszewskii (i2009i) . who suggested a dou- 
ble resonance 4:2:1 as the likely "survival recipe" for the entire 
three-planet system. 

As noted above, the major danger of system's destabilization 
comes from the two inner planets which are likely more massive 
and more tightly spaced. Keeping this in mind and taking into 
account that data on the innermost companion (e.g. its differen- 
tial proper motion) are as yet the least reliable, we checked the 
dynamical stability properties the system would have without 
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Fig. 13. Fraction of left-over planetesimals after 100 Myr 
of evolution in the outer {top) and inner {bottom) rings. 
Filled/hatched bars: nominal/low planetary mass. 

HR8799 d. We used the same setup for MERCURY6 and re- 
stricted our analysis to three cases: non-inclined configuration, 
/ - 15°, and / - 30°. The overall result is that, as expected, 
the absence of the inner companion would drastically improve 
stability: 

- The non-inclined configuration becomes stable over a period 
of 100 Myr — not only for the nominal masses, but also for 
much higher masses up to Mb - 22Mjup and Mc - 30Mjup. 
The rapid breakdown of the system in less than 10 kyr would 
only be guaranteed with masses as high as Mh = 33Mjup and 
M, = 45Mjup; 

- For i - 15° and all O - 0°-180°, the system with nominal 
masses is always stable over 100 Myr; 

- For i - 30°, the system with nominal masses is unstable for 
Q = 40°-75° and stable otherwise. 



5. Dynamical interaction between the planetesimal 
belts and the planets 

To check where the outer and inner planetesimal belts (main- 
taining dust rings d and b) would be truncated by planets, we 
have picked up one exemplary angle configuration from within 
the "stability spot" seen in Fig. [8] The point chosen is (/, Q) = 
(30°, 10°), which implies initial semimajor axes of 74.7, 41.8 
and 24.6 AU for HR8799 b, c and d, respectively. 

We launched 1000 massless planetesimals with uniformly 
distributed orbital elements {cj, = 0-0.2, ip = 0-10°, 
{LL>p,Q.p, Mp} = 0-360°). Of these, 200 planetesimals were ini- 
tially confined to a ring between Oinner = 10-20 AU and the 
other 800 between Oouter = 90 - 130AU. We integrated their 
orbits over 100 Myr with an accuracy of 10"'^ for low and nom- 
inal planetary masses, as given in Table |4] The high mass case 
was excluded, because the stability of companions themselves at 
(/, Q.) = (30°, 10°) was only marginal. 

Figure[T2]shows that both the inner part of the outer ring and 
the outer part of the inner ring are swiftly cleared by the adjacent 
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Fig. 12. Evolution of semimajor axes of planetesimals over 100 Myr of for the nominal {left) and low {right) planet masses. Starting 
planetesimal rings are shown with vertical bars. Solid and dashed lines represent the location of planets and important resonances, 
respectively. Note that the outer ring {top) is plotted in barycentric osculating elements and the inner ring in astrocentric ones. This 
choice leads to a sharper visibility of the resonance positions. 



planet (b and d, respectively). The belt of remaining planetesi- 
mals develops "Kirkwood gaps" at the positions of major reso- 
nances. Over tens of Myr, the gaps get progressively more pro- 
nounced. Simultaneously, the fraction of particles surviving be- 
tween them gradually decreases. As expected, the survival prob- 
ability of planetesimals at a given distance is larger for lower 
masses of the companions. 

For both rings and for the nominal- and low-mass cases. 
Fig. [13] presents the fraction of planetesimals that survived af- 
ter 100 Myr in orbits with different initial semimajor axes. For 
nominal planetary masses, the outer ring shows a considerable 
depletion, with only ~ 10 to 20% planetesimals surviving even 
outside 120 AU. There are almost no survivors inside the 5:3 and 
7:4 resonances with HR8799 b at 105-1 10 AU. In the low- 
mass case the survival fraction in the outer ring is appreciably 
higher (15-70% between k 1 10-130 AU). The inner ring retains 
80-100% of planetesimals inside 11 AU. The outer edge of the 
ring is at 13 AU (nominal-mass planet) to 15 AU (low-mass one). 

These results have to be compared with dust locations found 
from the SED fitting in section 3.2. As noted there, the far- 
infrared to millimeter part of the SED requires dust in the outer 
ring as close as 120 AU from the star; we have just shown that 
a significant fraction of planetesimals survives outside 120 AU 
after 100 Myr, at least in the low-mass case. In the nominal-mass 
case, the fraction of survivors is lower, but any firm conclusions 
appear premature, since the location of the outer ring is in fact 
not well-constrained (see section 3.2.). Next, the IRS spectrum 
interpretation requires dust in the inner ring at least at 10 AU 
away from the star This is comfortably within the stability zone 
inside the orbit of HR8799 d for both the nominal-mass and low- 
mass cases. Moreover, 10 AU quoted above is the distance where 



dust is required; as discussed in section 3.2, the parent planetes- 
imals would orbit closer to the star, being yet safer against the 
perturbations of the innermost planet than their dust. In sum- 
mary, our analysis of the outer system might slightly favor the 
low-mass case, but would not really pose any additional strong 
constraints to the planetary masses. 

6. Conclusions and discussion 

6.1. Conclusions 

In this paper, we made an attempt of a coherent analysis of var- 
ious portions of observational data currently available for the 
system of a nearby A5 star HR8799, which hosts debris dust 
as well as thre e planetary ca ndidates recently discovered via di- 
rect imaging dMarois et alj r2008). A dedicated analysis of all 
known components of the system (the central star, imaged com- 
panions, and dust) leads us to a view of a complex circumstel- 
lar system (Fig. [T4l i. It contains at least three planets in nearly- 
circular coplanar orbits bordered by two dust-producing plan- 
etesimal belts, one outside the planetary region and another in- 
side it. Each planetesimal belt is encompassed by a dust disk. 
The outer dust disk may have a considable extention, perhaps 
several hundreds of AU. 

Our specific conclusions are as follows: 

1. With previous estimates of stellar age ranging from 20 
to a; 1 100 Myr, high luminosity of observed cold dust may 
favor younger ages of ^ 50 Myr. A younger age would au- 
tomatically lower the masses for all three companions, es- 
timated through evolutionary models, making the masses 
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Fig. 14. A schematic view of the system HR8799. 



more consistent with dynamical stability results (see conclu- 
sion 5). 

2. The system is seen nearly pole-on. Our analysis of the stel- 
lar rotational velocity suggests an inclination of 13 - 30°, 
whereas i ^ 20° seems to be mandatory for the system to be 
dynamically stable (see conclusion 5). Thus we arrive at a 
probable inclination range of 20 - 30°. 

3. Our analysis of the available mid-infrared to millimeter pho- 
tometry and spectrophotometry data reveals the presence 
of two dust rings, and therefore two parent planetesimal 
belts, an "asteroid belt" at ~ 10 AU and a "Kuiper belt" at 
^ 1 00 AU. The dust masses are estimated to be » 1x1 0"^ 
and 4 x 10"^ for the inner and outer ring, respectively. 

4. Assuming that the system is indeed rather young ($ 50 Myr) 
and based on the ph otometry of the companions reported by 
iMarois et al.l (|2008), our estimates with several evolutionary 
models suggest the masses of the companions to be lower 
than 7 to lOMjup. 

5. We show that all three planets may be stable in the mass 
range suggested in the discovery paper by Marois et al. 2008 
(between 5 and 13Mjup), but only for some of all possible 
orientations. For (M/,, Mc, Mj) = (5,7,7)Mjup, an inclina- 
tion ; ^ 20° is required and the Une of nodes of the sys- 
tem's symmetry plane on the sky must lie within 0° to 50° 
from north eastward. For higher masses {Mb, Mc, M,/) from 
(7, 10, 10)Mjup to (11, 13, 13)Mjup, the consti-aints on both 
angles are even more stringent. The stability of the two inner 
planets is due to locking in the 2:1 mean-motion resonance, 
and the stability of the outer couple is supported by the 2:1 
commensurability, too. However, in many stable cases only 
one of the two resonances is strong. Another one is often 
shallow, with a circulating rather than librating resonant ar- 
gument. For "wrong" orientations, the stability only seems 
possible with planetary masses lower than most evolution- 
ary models would predict even for the youngest possible age 
(cf. Table |3]l. Should this be the case, this would necessitate 
revisions to the models. 



6. Both dust/planetesimal belts appear to be dynamically stable 
against planetary perturbations, provided the masses of com- 
panions are such that they themselves are dynamically stable 
against mutual perturbations. 

6.2. Prospects for future observations 

Given the paucity of observational data available to date, many 
of the quantitative estimates listed above are quite uncertain and 
should be taken with caution. However, there is little doubt that 
new observations will arrive soon, verifying these estimates and 
on any account reducing the uncertainties. 

Firstly, new observations of the planets themselves are ex- 
pected and would be of great value. Better astrometry, and there- 
fore a better determination of the orbits should become pos- 
sible with improving instruments and methods of astrometric 
observations, and because of the longer time spans. Next, new 
photometry observations are needed. The companions are de- 
tected so far in JHK and L'. The SED is relatively flat in this 
wavelength range for objects with temperature of roughly 500 
to 1500 K. Imaging photometric detections of the companions 
in either Gunn z (lyum) or M (5 fim) must be possible with 
8 to 10m class telescopes. They would allow one to constrain 
the objects' temperature, because the colors z-J and L-M and 
their differences depend on temperature strongly over the rele- 
vant range of 500 to 1500 K. Further, spatially resolved spec- 
troscopy of the companions may be possible with VLT/Sinfoni 
or Subaru/IRSC, but would be very challenging. If successful, it 
would place tighter constraints on temperature and gravity and, 
hence, radius and mass of the companions. 

Secondly, new data on the dust portion of the system would 
be particularly promising. For instance, a better mid-infrared 
photometry would result in more reliable dust mass and loca- 
tion of the inner dust belt ("exozodi"). One could also think 
of near- and mid-IR interferometry observations, which have 
proven very successful not only for exozodi studies, but also 
for stellar radius determination (see, e.g., JD i Folco et al. 200^ 
Abs il et al.ll2006HDi Folco e"tani2007l:lAbsiI et al.ll2008l) . While 
HR8799 is too faint to be observed with the presently oper- 
ating CHARA/FLUOR and Keck Interferometry Nuller instru- 
ments (e.g. for CHARA the K-magnitude of ^ 4 mag is needed, 
whereas HR8799 has K = 5.2 mag), this should become possi- 
ble in the near future, for instance with VLTI/PIONIER and the 
LBTI Nuller More observational effort is required for the outer 
disk as well. Resolving the outer debris disk, especially in scat- 
tered light, would answer several key questions at a time. On the 
one hand, it would further constrain the inclination of the entire 
system and the orientation of its line of nodes on the sky plane, 
drastically reducing the parameter space assumed in the dynam- 
ical simulations. On the other hand, the precise location of the 
inner rim of the disk could place a direct upper limit on the mass 
of HR8799 b, mu ch in the sa me way it was recently done for the 
Fomalhaut planet jKalas et al . 2008; Chiang et al. 2009). Since a 
high dustiness of the debris disk of HR8799 makes it a relatively 
easy target, whereas advantages of resolving it are obvious, suc- 
cess can be expected in the near futur^ 

Once the location and masses of the dust belts are better con- 
strained from observations, it will become possible to access the 
position, masses, and other properties of directly invisible plan- 



After submission of this paper, we became aware of recent suc- 
cessful observations: the outer disk of HR8799 was resolved at 24 and 
70 yum with Spitzer/MIPS (Kate Su, pers. comm). The data analysis in 
ongoing. 
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etesimal belts that produce and sustain that dust. This could be 
done with the help of elaborate collisional models (' Krivov et alj 
[2OOS) . The results could provide additional clues to the forma- 
tion history of the system. 



6.3. Origin and status of tine system 

The HR8799 system is among a few systems known to date 
to possess more than one planet and at l east one planetesi- 
mal/d ust belt. Other ex amples are HD3 8529 M oro-M artin et al. 
120071) and HD69830 (Beichma n et"ani2005HLisseet al.ll2007h . 
Further, it is not the only system with directly imaged compan- 
ions whose masses most likely fall into the "planetary" (as op- 
posed to brown dwarf) range. However, in some sense HR8799 
does appear unique for the moment. The orbits of companions 
extending up to x 70 AU are large, and their masses most likely 
too - probably almost at the limit of their dynamical stability 
against mutual perturbations. Even though we strongly argue 
that the masses are well below the deuterium burning limit, it 
is not clear whether the companions have formed in a "plan- 
etary" way (from the protoplanetary disk) or "stellar" way (as 
a multiple stellar system). In this sense, it remains questionable 
whether we are dealing with "true" planets. An argument in favor 
of the "stellar" formation would be, for instance, a low metal- 
Hcity (Fe/H ^ -0.47, Gray & Kaye 1999), atypical of — al- 
though not the lowest among — the known planet host stars. The 
low metallicity is particularly unusual for a system with several 
high-mass planets. "Planetary" way of formation in situ could 
be f easible through the gravitational instability (GI) (Cameron' 
11978 ; Boss 1998). For Gl to work, the densit y of the protop lan- 
etary disk should excess the Toomre density (lToomrelll964h . At 
the same time, it should be low enough to allow efficient cool- 
ing, which is required for a disk to fragment into bound clumps 
(lGammie 2001). Both radiative and convective cooling rates may 
not be efficient enough for direct formation of g iant planets by 
GI wi thin several tens of AU from the parent star (lRafikovll2005l 
l2007h . As far as the standard core accretion scenario is con- 
cerned, the only possibility to explain the formation of planets 
with several Jupiter masses would be to admit that they have 
formed closer to the st ar and then were scatter ed gravitation- 
ally to wider orbits (e.g. lVeras & Armitagell2004h . Nevertheless, 
it is difficult to find a mechanism that has circularized their or- 
bits subsequently. Alternatively, massive planets formed by core 
accretion could have smoothly migrated from their birth places 
outward. However, this would require displacing a comparable - 
and therefore an unrealistically large - mass of planetesimals 
inward over large distances, rendering this mechanism problem- 
atic. Thus, by and large it remains unclear how the planets have 
formed. The era of directly imaged extrasolar planets that has 
just begun must eventually bring answers to these and many 
other questions. 
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